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Reynolds-Averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) analysis 
was conducted to study the low-speed stall aerodynamics of a GUI aircraft’s swept wing 
modified with a laminar-flow wing glove. The stall aerodynamics of the gloved wing were 
analyzed and compared with the unmodified wing for the flight speed of 120 knots and 
altitude of 2300 ft above mean sea level (MSL). The Star-CCM+ polyhedral unstructured 
CFD code was first validated for wing stall predictions using the wing-body geometry from 
the First American Institute of Aeronautics and Astronautics (AIAA) CFD High-Lift 
Prediction Workshop. It was found that the Star-CCM+ CFD code can produce results that 
are within the scattering of other CFD codes considered at the workshop. In particular, the 
Star-CCM+ CFD code was able to predict wing stall for the AIAA wing-body geometry to 
within 1 degree of angle of attack as compared to benchmark wind-tunnel test data. Current 
results show that the addition of the laminar-flow wing glove causes the gloved wing to stall 
much earlier than the unmodified wing. Furthermore, the gloved wing has a different stall 
characteristic than the clean wing, with no sharp lift drop-off at stall for the gloved wing. 
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I. Introduction 

r J^HE National Aeronautics and Space Administration (NASA) Dryden Subsonic Research AircrafT (SCRAT) 
aircraft tail number 804 is pictured in Fig. 1. This modified Gulfstream Aerospace Corporation (Savannah, 
Georgia) GUI aircraft is being proposed for use as the flight testbed for the Subsonic Aircraft Roughness Glove 
Experiment (SARGE) [1] project. A collaborative effort between NASA and Texas A & M University (College 
Station, Texas), SARGE proposes to conduct flight research on the effectiveness of the discrete roughness elements 
(DREs) in increasing the laminar-flow region on a generic subsonic transport swept wing for flight conditions 
representative of modern subsonic jetliners. The SARGE project directly supports the goal of the NASA 
Environmentally Responsible Aviation (ERA) Project of enabling new technologies for the purposes of increasing 
aircraft efficiencies and reducing the adverse effects that aviation has on the environment. There are several reasons 
that SARGE is a good fit for the NASA Dryden SCRAT aircraft. First, the flight environment provides the very low 
freestream turbulence levels required to demonstrate the effectiveness of the DRE laminar-flow technology. Also, 
DREs have been found to produce significant increases in laminar-flow regions of a swept wing test article in flight 
at lower Reynolds numbers and Mach numbers [2], The Dryden SCRAT aircraft can help by expanding the flight 


research of the SARGE laminar-flow wing glove to higher Reynolds numbers and Mach numbers representative of 
typical full-scale jetliners. 




Fig. 1. NASA Dryden GUI SubsoniC Research AircrafT (SCRAT). 

conceptual sketch of the GUI aircraft modification is shown in Fig. 2. The SARGE DRE laminar-flow wing 
would be mounted only on the left wing of the SCRAT aircraft, resulting in an asymmetric aircraft 
configuration. The glove spans approximately 30 percent of the half span of the aircraft wing, occupying the 
midwing section, away from both the wing root as well as wing tip. The location of the glove was chosen to be away 
from the region of influence of the aircraft engine near the wing root for better flow quality over the glove, as well 
as away from the aircraft aileron, which is located towards the wing tip to avoid potential interference with the 
aerodynamic effectiveness of the left aircraft aileron. Micron-sized DREs mounted near the leading edge of the wing 
glove provide the spanwise-periodic roughness variation that controls the unstable crossflow vortices and delays the 
transition to turbulence on the wing glove. The actual DREs are not directly modeled in the present computational 
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fluid dynamics (CFD) simulations. 



Fig. 2. Conceptual sketch of the SARGE DRE laminar-flow wing glove mounted on the NASA Dryden 
SCRAT aircraft's left wing. 

All of the proposed SARGE flight-test points will be done around trimmed aircraft cruising flight conditions for 
the SCRAT aircraft. Hartshorn [3] has found that the wing glove produces very small changes in aircraft 
aerodynamics at aircraft cruising flight conditions. Laminar-flow wing glove flight testing has been conducted 
extensively in the past using a wide variety of aircraft such as the JetStar (Lockheed Martin, Bethesda, Maryland), 
Learjet (Learjet, now Bombardier Aerospace, Wichita, Kansas), Boeing 757 (The Boeing Company, Chicago, 
Illinois), F-lll (General Dynamics, now Lockheed Martin, Bethesda, Maryland), F-14 (Grumman, now 
Northrop-Grumman Corporation, West Falls Church, Virginia), F-15 (McDonnell Douglas, now The Boeing 
Company, Chicago, Illinois), and F-16XL (General Dynamics, now Lockheed Martin, Bethesda, Maryland). A 
preliminary literature search has found no mention of stalls as significant issues in all of these previous flight 
research efforts. The laminar-flow wing glove proposed for the left wing of the SCRAT does not represent a 
significant departure from these previous laminar-glove flight research efforts. Therefore, it is not expected that 
modified aircraft stall characteristics would present serious problems for conducting the SCRAT aircraft SARGE 
wing glove research flights at high-speed cruising flight conditions. However, low-speed stall aerodynamics of the 
gloved wing need to be analyzed in order to ensure flight safety for the takeoff and landing portions of the SARGE 
glove test flights. 

The laminar-flow wing glove introduces many significant physical modifications to the left wing of the SCRAT 


aircraft, as listed below: 


• Modifications of the basic wing airfoil section for 30 percent of the half span of the aircraft wing 
include: smaller leading edge (LE) radius, longer chord length, a smaller thickness over chord (T/C) 
ratio, and modified airfoil camber line. 

• There are two leading-edge snags introduced by the glove extending beyond the LE of the unmodified 
aircraft wing. 

• Four standard aircraft wing vortex generators (out of a total of 3 1) are removed to make room for glove 
installation on the left wing. The right wing retains the full 3 1 vortex generators. 

The DREs are expected to provide significantly more laminar flow over the wing glove than the equivalent area 
on the unmodified Gill aircraft wing. Also, the laminar-flow airfoil design of the wing glove modifies the transonic 
flow over the wing as compared to the normal GUI aircraft wing airfoil. For example, at an aircraft Mach number of 
0.75, the shock wave on the upper surface of the wing in the laminar glove area is further downstream compared to 
the shock location of the clean wing at the same flight Mach number. 

Finally, the glove is mounted only on the left wing, resulting in an asymmetric aircraft configuration. This 
asymmetric mounting of the glove on the left wing only could create undesirable asymmetrical forces and moments 
should the gloved left wing stall earlier than the unmodified right wing. The main objective of the current analysis 
effort is to determine if the addition of the laminar-flow wing glove to the left wing causes left wing stalls to occur 
before the unmodified right wing of the airplane. 

In this report, the Star-CCM+ polyhedral unstructured CFD code is first validated for low-speed wing stall 
predictions using the wing-body geometry from the First American Institute of Aeronautics and Astronautics 
(A1AA) CFD Fligh-Lift Prediction Workshop [4]. The effects of grid sizes and turbulence models on the accuracy of 
the CFD stall predictions are quantified. Finally, the low-speed stall aerodynamics of the gloved wing are presented 
and compared with the unmodified wing aerodynamics. 

II. CFD Code Validation Study 

Steady Reynolds- Averaged Navier-Stokes (RANS) CFD simulations have been found to be qualitatively 
accurate for predicting low-speed wing stalls in the First AIAA CFD High-Lift Prediction Workshop [4]. An 
assessment of the prediction capability of twenty-one participants from eight countries using a variety of different 


structured as well as unstructured CFD codes was made in this workshop. Although there was significant scattering 
in the CFD predictions as compared with the benchmark wind-tunnel test data, it was shown in [4] that the angle of 
attack (AoA) of the wing at stall could be predicted with reasonable accuracy. Since the CFD code Star-CCM+ was 
not considered in this workshop, it is necessary to conduct a CFD code validation study to ensure that the 
Star-CCM+ high-lift results are reliable. 

Figure 3 illustrates the First A1AA CFD High-Lift Prediction Workshop Configuration 1 wind-tunnel geometry 
that was used for the present validation study. In this figure, contours of pressure coefficient on the wing-body 
surface and the plane of symmetry are shown together with the flow streamlines, illustrating the post-stall separated 
flows over this A1AA wing-body configuration. The selected configuration is one of the two configurations 
considered by the A1AA workshop. In this configuration, the LE slat is deflected by 30 degrees and the flap is 
deflected by 25 degrees. The Star-CCM+ meshes were created incorporating the guidelines and best practices 
provided in [5,6]. Some of the more useful guidelines include: 

• Far-field located at a distance of -100 reference chords [5]. 

• Surface grid spacing should be approximately 0.3 percent of the airfoil chord [6]. 

• Initial y+ wall spacing normal to the wing surface should be less than one [5] or equal to one [6]. 

• Geometric stretching ratio not more than 1.2 [6]. 



Fig. 3. Geometry for the First AIAA CFD High Lift Prediction Workshop configuration 1. 


The Star-CCM+ code validation study was conducted at a free-stream Mach number of 0.2 and a Reynolds 
number based on the mean aerodynamic chord of 4.3 million as specified by the A1AA workshop. All runs were 
conducted frilly turbulent. From a turbulence models sensitivity study, it was found that, using the same mesh with 
the same initial y+ normal wall spacing values less than 1, the automatic y+ turbulence modeling option in 
Star-CCM+ caused the wing-body CFD simulations to stall too early, and the low-y+ option produced more accurate 
stalled AoA predictions in better agreement with the A1AA workshop results. 

As can be seen in Figs. 4 to 6 below; the lift, drag, and pitching moment coefficients computed from the 
Star-CCM+ CFD simulations quantitatively agree with the wind-tunnel data from the First AIAA CFD High-Lift 
Workshop. Results using the Spalart-Allmaras turbulence model are shown for three different mesh sizes: coarse 
(22 million cells), medium (30 million cells), and fine (46 million cells). It can be seen that the Star-CCM+ results 
are within the bounds of scattering of the results from most of the CFD codes used in this workshop. There was even 
more scattering in the CFD codes considered in the workshop comparison than indicated by the green boundaries in 
Figs. 4 to 6, and the green boundaries in these figures exclude the outliers found in the workshop’s results. 

■ Wind tunnel data 
— © — Star-CCM+ SA coarse (22e6 cells) 

— A — SA medium (30e6 cells) 

— I — SA fine (46e6 cells) 


AIAA Hi-Lift CFD codes scattering 



Fig. 4. Star-CCM+ CFD code validation results for lift coefficient, Mach 0.2, Reynolds number based on 
mean aerodynamic chord of 4.3 million. 
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Fig. 5. Star-CCM+ CFD code validation results for drag coefficient, Mach 0.2, Reynolds number based on 
mean aerodynamic chord of 4.3 million. 
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Fig. 6. Star-CCM+ CFD code validation results for pitching moment coefficient, Mach 0.2, Reynolds number 
based on mean aerodynamic chord of 4.3 million. 

Although there are disagreements with the wind-tunnel test data, it can be seen that the Star-CCM+ CFD results 
are within the scatterings of most of the CFD codes considered in the workshop and tend to follow the trends of the 
other CFD codes well. For example, most of the codes including Star-CCM+ under predict both the drag coefficient 



as well as the pitching moment coefficient (less negative), as shown in Figs. 5 and 6, respectively. Most importantly, 
for the objective of the current work, it can be seen in Fig. 4 that Star-CCM+ is capable of predicting wing stall to 
within 1 degree of AoA of the wind-tunnel test data. The wind-tunnel lift coefficient data in Fig. 4 shows that wing 
stall occurs at an AoA of 36 degrees. The coarse-grid Star-CCM+ CFD result shows a stall AoA value of 
35 degrees. Both the medium- and fine-grid Star-CCM+ results show a stall AoA of 37 degrees, indicating grid 
convergence for this particular result. 

III. Laminar-Flow Gloved Wing Stall Analysis 

Figures 7(a) and (b) show the CAD models for the Star-CCM+ CFD swept wing simulations. The clean wing is 
shown in Fig. 7(a), and the gloved wing is shown in Fig. 7(b). Only one wing on a plane of symmetry was modeled 
in the current CFD simulations to concentrate as many grid cells as possible on the wing. The vortex generators can 
be seen as a row of small dots toward the wing tip in these figures. For clarity, only the mesh on the symmetry plane 
is displayed in Figs. 7(a) and (b). There is such fine mesh spacing on the wing surfaces that these wings would 
appear solid black if the wing mesh is turned on in Figs. 7(a) and (b). Figures 8(a) and (b) illustrate a closed-up view 
of both wing configurations in the vicinity of the vortex generators near the wing glove location, with the surface 
mesh turned on. It can be seen that four vortex generators are missing in the gloved wing in Fig. 8(b) as compared to 
the clean wing in Fig. 8(a). The appropriate number of vortex generators was gridded for both the clean wing 
(31 vortex generators total) as well as the gloved wing (27 vortex generators total). As discussed in the introduction 
section above, four standard aircraft wing vortex generators are removed to make room for glove installation on the 
left wing. Two different mesh sizes were considered: medium (35 million cells) and fine (45 million cells). All of 
these results were obtained using the low-y+ turbulence modeling option in Star-CCM+ with the first y+ value from 
the wall being approximately 0.1. It was found in the CFD code validation study section above that this y+ option 
produces the most accurate stall results. All CFD runs were conducted fully turbulent. Although there is flow 
transition on the laminar glove, it is likely that the flow rapidly transitions to mostly turbulent as soon as the AoA is 
increased above the design angle for the laminar-flow glove due to the off-design pressure gradient over the upper 
surface of the glove at higher AoA. 

Approximately 19 prism layers were used in the boundary layer over the wing with an overall prism layer 
thickness of 0.055 m. Far-field boundaries are placed approximately 100 times the mean aerodynamic chord (MAC) 


away from the wing. For the medium mesh, surface target mesh size is 0.3 percent MAC with a minimum surface 
mesh size of 0.075 percent MAC. Both of the surface mesh grow rate and volume cell grow rate are set at 1.2. All of 
the simulations were started at a low AoA of just below 6 degrees then gradually increasing until past lift stall. The 
converged solution at a lower AoA is used to start the next higher AoA calculation. Convergence is achieved when 
forces and moments, as well as all of the residuals become either constant for fully attached wing flows or periodic 
for the separated, post-stalled wing flows. 

Figure 9 summarizes the Star-CCM+ CFD wing lift coefficient stall results for the clean wing as well as the 
gloved wing for the low-speed, low-altitude flight conditions of 120 knots and 2300 ft above mean sea level (MSL), 
which is approximately ground-level at Edwards Air Force Base (Edwards, California). Before stalls, the gloved 
wing is seen to have slightly higher lift than the clean wing. Hartshorn [3] has shown that the airplane can 
compensate for this small difference in lift with no problems. The gloved wing stalls much earlier than the clean 
wing. Furthermore, the stall characteristic of the gloved wing is different from the clean wing. When the clean wing 
stalls at an AoA of approximately 17 degrees, there is an abrupt loss of lift. However, when the gloved wing stalls at 
an AoA of approximately 10.5 degrees, there is no abrupt loss of lift. Instead, the lift remains more or less constant 
at its value at AoA of 10.5 degrees. Although not captured in these steady RANS CFD simulations, the post-stall 
flows over the wing are likely very unsteady for both the clean as well as gloved wings. 



a) Clean wing. 


Fig. 7. Star-CCM+ CFD CAD geometry for the clean wing and gloved wing. 



b) Gloved wing. 
Fig. 7. Concluded. 




b) Gloved wing (4 vortex generators missing). 

Fig. 8. Closed-up view of the Star-CCM+ CFD grid for the clean wing and gloved wing with vortex 
generators. 
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O-deg flap gloved wing (medium CFD grid) 

• Gloved wing (fine CFD grid) 



Fig. 9. Star-CCM+ CFD wing lift coefficient stall results, 120 knots, 2300 ft above MSL. 

Four vertical gray dotted lines in this figure denote the trimmed AoA for four flap configurations for the normal 
(ungloved) wing of the clean airplane. When no aircraft flap is used (0-degree flap), it can be seen that the trimmed 
AoA of the clean airplane is actually higher than the stall AoA of the gloved wing. Since the left gloved wing stall 
occurs before the normal aircraft trimmed AoA, asymmetrical wing stall can occur without warning if no flap is 
used. When the aircraft flaps are deflected to 10, 20, or 39 degrees, the trimmed AoA values of the aircraft can be 
reduced to below the gloved wing stall AoA, possibly allowing aircraft operation at this low-speed, low-altitude 
flight condition. Note that the laminar-flow glove would likely change the stall characteristics of the wing with flap 
deflected, and an analysis of gloved wing stalls with flap deflections will be needed before flap deflections can be 
considered as mitigation. In addition to flap deflections, higher takeoff and landing speeds could be considered to 
increase the gloved wing stall margins. 

Figure 10 compares the flows over the clean wing, Fig. 10(a), and the gloved wing, Fig. 10(b), at the aircraft 
AoA of just below 12 degrees. As shown in Fig. 9 above, at this AoA the clean wing still provides lift with plenty of 
stall margin while the gloved wing has stalled. It can be seen that the flow remains smooth and fully attached for the 
clean wing, while for the gloved wing; the flow over the glove is separated causing the gloved wing to stall earlier 
than the clean wing. A possible explanation for the different stall characteristic of the gloved wing seen in Fig. 9 is 
provided by examining Fig. 10(b). Only the flow over the glove and its vicinity is stalled, with the rest of the flows 


over the wing remaining attached and providing some lift. Therefore, there is no sharp lift drop-off at gloved wing 
stall as compared to the clean wing stall. 
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a) Unstalled clean wing. 
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b) Stalled gloved wing. 

Fig. 10. Comparison of flows over unstalled clean wing versus stalled gloved wing at an AoA of just below 12 
degrees, 120 knots, 2300 ft above MSL. 

IV. Conclusions and Recommendations 

In conclusion, Reynolds-Averaged Navier-Stokes computational fluid dynamics (CFD) analysis was conducted 


to study the low-speed stall aerodynamics of a swept wing modified with a laminar-flow wing glove. The stall 



aerodynamics of the gloved wing were analyzed and compared with the unmodified wing at one low-speed and low- 
altitude flight condition. First, the Star-CCM+ polyhedral unstructured CFD code was validated for wing stall 
predictions using the wing-body geometry from the First A1AA CFD High-Lift Prediction Workshop. It was found 
that the Star-CCM+ CFD code can produce results that are within the spread of other CFD codes considered at the 
workshop. In particular, the Star-CCM+ CFD code was able to predict wing stall for the A1AA wing-body geometry 
to within 1 degree of angle of attack (AoA) as compared to the benchmark wind-tunnel test data. Current analysis 
results show that the addition of the laminar-flow wing glove causes the modified and gloved aircraft wing to stall 
much earlier than the unmodified wing. The gloved wing also has a different stall characteristic than the clean wing, 
with no sharp lift drop-off at stall. These results have important flight safety consequences, since the left gloved 
wing could stall prematurely before the right clean wing, producing a serious asymmetric aircraft configuration 
issue. It is possible to use flap deflections in order to reduce the aircraft’s trimmed AoA to values below the stall 
AoA of the gloved wing, potentially allowing the laminar-flow gloved wing aircraft to operate at the low-speed 
flight condition under consideration. However, an analysis on gloved wing stalls with flap deflections needs to be 
done before this can be considered as mitigation. In addition, higher aircraft takeoff and landing speeds might be 
used to increase the gloved wing stall margin. 

Only 0-degree flap was considered in this work. The addition of the laminar-flow wing glove could change the 
stall characteristics of the wing with deflected flaps. Also, ground effects were not considered in this study and 
could potentially modify the stall characteristics of the gloved wing. Both of these issues should be analyzed in 
future studies before the gloved wing aircraft could be flown. 

It is further recommended that the modified gloved wing lift, drag, and moment results obtained in this report 
need to be used in a 6-degrees of freedom GUI flight simulation software to determine gloved wing aircraft stall 
margins during takeoffs and landings. 
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